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We present a theoretical study of in-plane magnetization 
reversal for vicinal ultrathin films using a one-dimensional mi- 
cromagnetic model with nearest-neighbor exchange, four-fold 
anisotropy at all sites, and two-fold anisotropy at step edges. 
A detailed "phase diagram" is presented that catalogs the 
possible shapes of hysteresis loops and reversal mechanisms 
as a function of step anisotropy strength and vicinal terrace 
length. The steps generically nucleate magnetization reversal 
and pin the motion of domain walls. No sharp transition sepa- 
rates the cases of reversal by coherent rotation and reversal by 
depinning of a 90° domain wall from the steps. Comparison 
to experiment is made when appropriate. 

PACS numbers:75.70.Ak,75.60.-d,75.10.Hk 



I. INTRODUCTION AND BACKGROUND 

Laboratory studies of ultrathin films of transition met- 
als confirm the general principle that broken symmetry 
induces magnetic anisotropy. |l],|^] The most common ex- 
ample is the loss of translational invariance at the free 
surface of a film or at an internal interface of a multilayer 
structure. The phenomenological "broken-bond" model 
of Neel Q then provides an intuitive way to understand 
why atoms at the surface or interface favor alignment of 
their magnetic moments either parallel or perpendicular 
to the broken symmetry plane. [Q In some cases, perpen- 
dicular anisotropy occurs that is strong enough to over- 
whelm the tendency for in-plane magnetization favored 
by magnetostatic shape anisotropy. This situation can 
be exploited for a variety of applications and has been 
the subject of very thorough experimental and theoreti- 
cal work. H 

In this paper, we focus on a related phenomenon: the 
magnetic anisotropy induced by crystallographic steps on 
the surface of a single crystal film. Here, it is the loss 
of translational invariance in directions parallel to the 
(nominal) surface plane that is germane. Application of 
the Neel model suggests that local moments will tend to 
align themselves either parallel or perpendicular to the 
local step orientation. The magnitude of the effect (on a 
per atom basis) is predicted to be comparable to conven- 
tional surface anisotropy. However, it was not until 1987 
that Hillebrands et al. invoked step-induced anisotropy 
to rationalize their surface spin wave data for epitaxial 
Fe/W(110). [§ Since all ultrathin films invariably have 
step edges (associated either with steps on the substrate 
or with the nucleation and growth of monolayer height 



islands during the growth process) it is not surprising 
that subsequent experimental studies often cite this phe- 
nomenon in connection with "surface roughness effects" . 
I 

We recently presented a theoretical study of in-plane 
magnetization reversal in ultrathin films with step struc- 
ture in typical samples. Q The model film was com- 
prised of an array of square, monolayer-height, magnetic 
islands of variable size and density on top of a few com- 
plete magnetic layers. Classical XY-type spins at each 
site were presumed to rotate in the surface plane subject 
to nearest-neighbor ferromagnetic exchange, an intrinsic 
four-fold in-plane anisotropy at all surface sites, Zeeman 
energy from an external field, and a two-fold anisotropy 
at island perimeter sites only. Numerical simulations and 
simple geometric scaling arguments predicted significant 
variations in coercivity as a function of coverage for layer- 
by-layer growth at low island nucleation densities. This 
result was found to be in semi-quantitative agreement 
with the surface magneto-optic Kerr effect (SMOKE) 
data of Buckley et al. |] for the Cu/Co/Cu(001) sys- 
tem. A subsequent Monte Carlo simulation study 
of coercivity in islanded Fe sesquilayers on W(110) using 
an in-plane Ising-type spin model yielded similarly good 
results in comparison to experiment. 

The theoretical results of Ref. || were interpretable on 
the basis of several qualitative concepts: (i) nucleation of 
magnetization reversal at island edges; (ii) pinning of do- 
main walls at island edges; and (iii) fusion of nearby do- 
mains. Unfortunately, even the simple island morphology 
studied there was still too complex to permit a detailed 
analytic treatment of the reversal process as one might 
desire. For this reason, we analyze an even simpler prob- 
lem in this paper: zero-temperature, in-plane magnetiza- 
tion reversal in ultrathin vicinal films. The basic model 
sketched above remains unchanged except that the mor- 
phology is simplified to a periodic array of flat magnetic 
terraces separated by straight, monolayer-height steps. 
This renders the problem one-dimensional and amenable 
to analytic study. 

One-dimensional models of magnetization reversal 
with inhomogeneous or competing anisotropics have been 
a fixture of the magnetism literature for many years. 
Most of these papers focus on the demonstration that 
planar defects in bulk ferromagnets can nucleate rever- 
sal and/or pin domain wall motion. If operative, these 
effects call into question the suitability of the popular 
single-domain, coherent rotation model of Stoner and 
Wohlfarth as a description of magnetization rever- 
sal. Filipov nff2| and later Brown p3| studied the effect 
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of surface anisotropies on the nucleation field (where the 
magnetization first deviates from its saturation value) 
while Mitsek and Semyannikov |TJ| and later Friedberg 
and Paul [15| focused on the depinning of pre-existing 
reversed domains as a determinant of the coercive field 
(where the magnetization projected on the external field 
direction first falls to zero). In recent years, Arrott has 
been explicit in the application of these ideas to ultrathin 
films with and without step structure. p6[ Our analysis 
will be seen to substantially extend all of these studies. 

On the experimental side, Heinrich et al. [[l7| first 
drew attention to the fact that a step-induced uniaxial 
anisotropy must be present on vicinal surfaces. Subse- 
quent work confirmed this observation ||l^-pO|] and re- 
vealed a number of other systematic features. As particu- 
lar motivation for the present work, we draw attention to 
the SMOKE data of Kawakami et al. obtained from 
Fe films grown on stepped Ag(OOl) substrates. Charac- 
teristic "split-loop" hysteresis curves were found where 
the degree of splitting varied smoothly with the degree 
of vicinality. The authors interpreted their results us- 
ing a single domain switching model where the step edge 
anisotropy was distributed over the entire surface. The 
analysis below will make clear the extent to which this 
description can be regarded as reliable. 

The plan of our paper is as follows. Section II is an 
overview that includes (i) a discussion of the model as- 
sumptions; (ii) the definition of important dimensionless 
quantities and the presentation of a "phase diagram" that 
catalogs the possible hysteresis loop topologies than can 
occur; (iii) a qualitative discussion of the physical mech- 
anisms of magnetization reversal that can occur; and (iv) 
a preliminary comparison to relevant experiments. Sec- 
tion III reports our mathematical procedures. We define 
the Hamiltonian used and solve the model exactly to ex- 
tract the physics of zero-temperature reversal in the sin- 
gle domain and single step limits. The intermediate case 
of multiple steps is formulated and solved numerically. 
Section IV is a discussion that complements the earlier 
overview in light of our analytic and numerical results. 
We consider the crossover between coherent rotation and 
domain wall depinning, discuss relevant experiments in 
more detail, and comment on various limitations and ex- 
tensions of the model. Section V summarizes our results 
and concludes the paper. 



II. OVERVIEW 
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FIG. 1. Geometry and anisotropies for a monolayer of mag- 
netic material on a vicinal non-magnetic substrate. The sub- 
strate steps are periodically separated by a distance L. There 
is a four-fold anisotropy everywhere on the surface, and a 
strong two-fold anisotropy localized at the steps. 



flat terraces of length L separated by monoatomic height 
steps. We assume perfectly straight steps so that the 
spin configuration is a function only of the spatial coordi- 
nate perpendicular to the steps. The problem is thereby 
reduced to a one-dimensional classical spin chain with 
ferromagnetic exchange J. 

The total surface anisotropy from all sources is pre- 
sumed to compel the spins to lie in the plane of the 
substrate terraces. To model surfaces with cubic sym- 
metry, we assign a four-fold anisotropy with strength K± 
to every site of the chain and a two-fold anisotropy with 
strength K 2 to every step site. The sign of K4 is chosen 
to favor spin orientations parallel and perpendicular to 
the steps. 23 If the sign of K 2 favors spin orientation 
parallel (perpendicular) the steps, we apply the external 
field H perpendicular (parallel) to the steps. These cases 
are identical by symmetry. Magnetostatics contributes 
to the total surface anisotropy that compels the spins to 
lie in-plane. For this model, with in-plane spins, magne- 
tostatics is not treated explicitly because its additional 
effects arc known to be negligible in the ultrathin limit. 
p4j Figure |l| is a schematic representation of the physical 
situation and the spin chain model studied here. 

For simplicity, we choose units where the lattice con- 
stant a is one and J, K 2 , and H, all have units of 
energy. To recover dimensional units as used in B, di- 
vide K4 and K 2 by a 2 , and divide H by /i, where fi is 
the atomic magnetic moment. 



A. Model Assumptions 



We consider a uniformly thick ultrathin magnetic film 
adsorbed onto a vicinal non-magnetic substrate. By flat, 
we mean that the film has no island structure, e.g. a 



film grown in step flow mode. |22| By ultrathin, we mean 
that there is no significant variation in the magnetization 
density in the direction perpendicular to the plane of the 
substrate terraces. By vicinal, we mean a sequence of 



B. The Phase Diagram 

We organize our discussion of hysteresis in this system 
around a "phase" diagram (Figure |^) whose axes are a 
scaled step anisotropy strength JC = K 2 /2a and a scaled 
terrace length C = L/W where K 2 is the step anisotropy 
energy, a = \J2JK± is the domain wall energy, and W = 
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\J J/2K4 is the exchange length. The solid lines delineate 
four distinct hysteresis loop topologies. The dashed lines 
divide Phase II into three sub- variants. 

Figure [3] illustrates representative hysteresis loops in 
each phase. Since all the loops are symmetric with re- 
spect to the sign of H, it will be convenient to restrict 
discussion to the situation where the field changes from 
positive to negative. We define three characteristic values 
of the external field. The first deviation of the magneti- 
zation from saturation occurs at the nucleation field H^ ■ 
A jump in magnetization that initiates at the steps is 
denoted Hs- A magnetization jump that initiates on the 
terraces is denoted —Ht- Hs = Hm in Phases lie, III & 
IV. 

In phase I, all spins rotate continuously from the satu- 
ration direction to the reversed direction as the exter- 
nal magnetic field is reversed adiabatically. Near the 
left hand side of the phase diagram, the spins rotate 
nearly coherently as a single unit. This is called Stoner- 
Wohlfarth behavior. jll| But near the right hand bound- 
ary of the Phase I field, the spins near the step edge 
rotate more (per unit change in external field) than do 
the spins near the center of the terrace. There is no hys- 
teresis, i.e., no jumps appear in the magnetization curve, 
merely more or less spatially inhomogeneous spin rota- 
tion. 

In Phase Ila, spins within an exchange length of a step 
rotate away from the saturation direction at Hjy in re- 
sponse to the torque applied by the step anisotropy. A 
domain wall thus forms between the step spins and the re- 
maining terrace spins. A field-dependent energy barrier 
Aow separates this configuration from a configuration 
where all spins point nearly 90° from the saturation di- 
rection. Ajjw — > at H$ and the domain walls "depin" 
from the steps and sweep across the terraces. The accom- 
panying jump in magnetization is followed by a continu- 
ous segment of the hysteresis curve that passes through 
the origin. This is an SW-like regime of nearly coherent 
spin rotation. During this rotation, an energy barrier 
Asw separates the terrace spin configuration from the 
nearly reversed state. At H = —Ht, Asw disappears 
for the terrace spins farthest from the steps and a sec- 
ond jump in magnetization occurs. Reversal completes at 
—Hn when the step spins finally complete their rotation. 

Phase lib differs from Phase Ila because Ht > -Hjv 
and the final jump in magnetization carries the system di- 
rectly to the saturated reversed state. The phase bound- 
ary is the locus of points where Ht — H^. Note that 
there is a small range of K, where one encounters the 
phase sequence Ila — > lib — > Ila as C decreases from 
large values. 

Phase lie mostly occupies a portion of the phase di- 
agram where ICC < 1. Is this regime, the independent 
domain wall description used above is no longer appropri- 
ate because the walls have overlapped to the point where 
the magnetization inhomogeneity across each terrace is 
not large. The reversal is better described as nearly co- 
herent rotation, as above, where the degree of rotation 
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FIG. 2. Loop structure phase diagram. The independent 
variables are a scaled two-fold anisotropy strength at the step, 
K. and a scaled step separation C. Roman numerals label four 
distinct loop topologies. Lower case letters label three vari- 
ants of Phase II. The vertical and horizontal arrows respec- 
tively show the JC — > 00 and C —* 00 limits of the nearby 
phase boundaries. 

differs for spins near and far from the steps. On the 
other hand, a thin sliver of the lie phase field extends 
to very large values of C where the independent domain 
wall picture remains valid. This shows that there is no 
rigid correspondence between phases and reversal mech- 
anisms. More typically, as in this case, there is a smooth 
crossover from a domain wall picture to a coherent rota- 
tion picture. 

Phase III occupies the smallest portion of the phase di- 
agram. The step anisotropy here is sufficiently small that 
a negative field is needed to nucleate reversal. Otherwise, 
the reversal mechanism is identical to Phase lie. 

Phase IV is characterized by Hm < —Ht so that only 
a single magnetization jump occurs. In fact, Hn is so 
negative that the state with terrace spins nearly parallel 
to the step is not stable as it was in Phase III. During 
the jump, the degree of spatial homogeneity of the spin 
rotation is dictated by the magnitude of ICC. Nearly 
coherent SW reversal occurs when ICC <gC 1 while rotation 
initiates at the step when ICC >• 1. 

C. Relevant Experiments 

Two recent experimental studies of magnetization re- 
versal in thin iron films deposited onto vicinal and 
(nominally) flat surfaces can be interpreted with our 
phase diagram. Chen and Erskine JlSj ] studied ultrathin 
Fe/W(001) where the step anisotropy favors magnetic 
moment alignment perpendicular to the step. Their re- 
sults for an external magnetic field aligned parallel to 
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FIG. 3. Hysteresis loops. The type of hysteresis loop in 
different parts of the phase diagram, (see Fig. The scaled 
parameters for each loop are: I, fC = 1.25, £ = 0.5; Ha, 
K = 1.25, £ = 2.0; lib, K = 1.1, £ = 2.0; lie, K = 0.5, 
£ = 0.25; III, K = 0.5, £ = 0.75; IV, K. = 0.5, £ = 2.0. 



the steps can be compared with our results by symmetry. 
They observe loops characteristic of Phase III and Phase 
II for the samples they label "smooth" and "stepped" for 
1.5 ML iron coverage. 

Kawakami et al. pjj presented a sequence of four hys- 
teresis loops for the Fe/Ag(001) system that we interpret 
similarly as a transition from Phase III to Phase II. In 
this case, the step anisotropy favors magnetic moment 
alignment parallel to the step and the data they present 
for the external field aligned perpendicular to the step is 
relevant. More details of this comparison can be found 
in the Discussion section. 



III. QUANTITATIVE ANALYSIS 

A. General Results 

In the continuum limit, the model assumptions stated 
at the beginning of section II lead us to the following 
expression for the magnetic energy per unit length of 
step for an ultrathin film on a vicinal surface: 



E 



dx 



1 (d0 

2 \dx 



-Kacos49 — Hcost 



-K2 S(x — xs) cos 29 



(1) 



We remind the reader that J, K4, K2, and H all have 
units of energy. The lattice constant is unity so the inte- 
gration variable x is dimensionless. The function 9(x) is 
the angular deviation of the magnetization density from 
the field direction at point x. For definiteness, we take 
the latter to be perpendicular to the steps and pointing 
down the vicinal staircase of Figure [l| Note that the 
two-fold anisotropy acts only on step edge spins at the 
discrete positions xs- 

We seek spin configurations 6(x) that correspond to 
local minima of ([!]). In general, an energy minimum 
moves smoothly in configuration space as H changes and 
the corresponding spin configuration and magnetization 
change smoothly as well. Apart from accidental degen- 
eracies, the only exception to this behavior occurs when 
the energy minimum evolves to a saddle point. At that 
point, the spin configuration changes discontinuously, a 
new energy minimum is adopted, and a jump appears 
in the magnetization curve. Our goal is to calculate the 
field values where these jumps occur. Their number and 
sign distinguish the phases of the system. 

The Euler-Lagrange equation that determines the ex- 
tremal configurations of (|l|) is pq] 

d 2 9 

J— r = H sin 9 + 2K 4 sin 40 - V 5(x - x s )K 2 sin 29. 



dx 



s 



(2) 
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We seek solutions of this equation with the same peri- 
odicity as the steps. These solutions are parameterized 
by two constants, the spin angle at the center of each 
terrace 9t and the spin angle at each step 9$. One equa- 
tion that relates these two is obtained as follows. Place 
the origin x — at a step, multiply ^ by d9/dx, and 
integrate from the center of the terrace (x = —L/2) to 
an arbitrary point x on the same terrace. The result is 

1 1 / d9 \ 2 
Hcos9 T + -Ka cos46» T = -J — 

2 2 \dxJ 

+ Hcos9 + ~K 4 cos 46* (3) 

using the fact that d9 jdx = at the center of the terrace. 

The constant 9s appears when we evaluate (^) at a 
step. For this purpose, integrate (^) from x = 0~ to 
x = + and use reflection symmetry across the step, i.e., 



to get 



dx 



dx 



o+ 



d6 

dx 



K 2 sin26> 5 . 



(4) 



(5) 



Substitution into 



K 2 sin 2 26» s 



o- 

yields 

H(cos9t — cos#s) 

+ (cos4# T - cos 49s)/ 2 



(6) 



which relates 9t and 9s as desired. The scaled magnetic 
field H = H/K A . 

A second relation between 9t and 6*5 can be found 
that involves the terrace length explicitly by integrating 
(||) from the center of a terrace to the step edge: 



C = 2 



9 S 



d9 



e T y/[H(cos6 T - cos 9) + (cos46» T - cos40)/2] 



(7) 



The analysis to this point is completely general and 
forms the basis for all the approximate analytic and nu- 
merical results that follow. We begin our discussion with 
two special situations that can be treated in full analyt- 
ically: the single domain limit and the single step limit. 



where the value of the effective two-fold anisotropy K 2 — 
K 2 /L, as can be verified by substitution of a uniform 
spin configuration 9{x) = 9 into ([!]). 

In terms of the magnetization M = cos 9, we seek the 
stationary points of the quartic expression 

E = -K 4 (2M 2 - l) 2 + K 2 M 2 - HM, (9) 

i.e, the solutions of 



dF 1 

— =sm9[H- H(M)\ = 0. (10) 
d9 



where 



H(M) = (2K 2 + 8K 4 )M - 16K 4 M 3 . (11) 



The extremal condition is satisfied trivially when the 
magnetization is parallel or antiparallel to the field di- 
rection where sin6> = 0. But it is also satisfied by the 
cubic equation H = H(M). In either case, we must have 



d 2 E 



cos 9[H - H{M)\ + sin 2 9 



dH(M) 

d~M 



> (12) 



to guarantee that the solution is a local minimum of the 
energy. 

The first term on the right hand side of ( |l2"| ) determines 
the extremal properties of the sin# = solutions. The 



solution is a local minimum for H > where 



Hg = 2K 2 



8K 4 



(13) 



is the limiting value of the nucleation field when CJC <C 1 . 
Notice that portions of Phases I, lie, III, and IV appear 
in this limit where the nucleation field and the first 
jump field Hg are coincident. The 9 = ir solution is a 
local minimum for H < —H^. 

The second term on the right hand side of (|l^) deter- 
mines the extremal properties of the H = H(M) solu- 
tions. Because the coefficient of the cubic term is nega- 
tive, at most one of the three solutions to the cubic equa- 
tion satisfies dH(M)/dM > 0. This means that the mag- 
netization increases (decreases) when the field increases 
(decreases)-a condition that is met when \H\ < \Hj.\ 
where 



K2 
AKt 



,)5/2 



(14) 



B. The Single Domain limit 

This section focuses on the bottom left corner of the 
phase diagram where CIC -C 1. This is the Stoner- 
Wohlfarth limit where only a single homogeneous mag- 
netic domain is present. The energy per terrace per unit 
length of step E = E/L is 

E = -^K 4 cos 49 +^K 2 cos29-H cos 9, (8) 



is the limiting value of the jump field when CIC <C 1. 
This is true unless K 2 > 20K 4 in which case the H = H 
solution is stable for all values of M and there are no mag- 
netization jumps for any value of external field. When 
K 2 + 8K 4 < 0, the H = H{M) solution is never stable 
and the sin 9 = solutions are the only local minima. 

The above results can be applied to find analytic for- 
mulae for the three phase boundaries in the lower left cor- 
ner of the phase diagram. The system is in phase I when 
K 2 > 20K 4 since, as noted, the magnetization curve has 
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no jumps. The remanent slope is dM/dH = \j{2K% + 
8-K4). For 20i^4 > K2 > the system is in phase 

He. The remanent slope is dM/dH = l/(2K 2 + 8K 4 ). 
For 4K4 > K2 > 2K4 the system is in phase III. Phase 
IV occurs when 2K4 > Using these results and 

K2 = K2/L, the boundaries between the phases near 
the origin are: JC — 5C between phases I and lie, tC = £ 
between phases He and III, and JC = \L between phases 
III and IV. 



C. The Single Step limit 

The right edge of the phase diagram where C — > 00 is 
the limit where the step separation is large compared to 
the exchange length and the (somewhat larger) domain 
wall width. In that case, it is sufficient to study the case 
of a single step bounded by semi-infinite terraces on each 
side. Our goal again is to calculate the nucleation field 
i?7V and the jump fields H$ and Ht- We do this by 
focusing attention on the spin at the step where 9 = 8s 
and the spins at ±00 where we assume that 9 approaches 
the constant value 9t- 

The fact that 9(x) — > 9t as x — > ±00 implies that all 
spatial derivatives of 6(x) vanish at infinity. Applying 
this to (0) yields 



7r/2 and M = 0. An explicit formula for H^? can be 
found by noting that this jump initiates with the terraces 
spins at ±00. These obey the pure Stoner-Wohlfarth 
dynamics of Section II. B with K2 = 0. In particular, 
( |l5| ) is identical to (|l0|). The final magnetization jump 
thus occurs at 



-Hj? where 



m 



8V6 



Ka. 



(20) 



This value is a lower bound for the jump field when the 
terrace length is finite because the presence of nearby 
steps retards the final transition to the reversed state. 

The case < applies to Phases III and IV. The 
above discussion shows that at nucleation in Phase III, 
the saturated state jumps immediately to the spin con- 
figuration that satisfies ( [l5|) with 9t ^ 0. This state 
evolves smoothly until the magnetization jump at —H^?. 
In Phase IV, there is only a single jump because now 
( |i"5| ) has stable solutions only at 9 = and 9 — tt when 
H 



TTOO 

N ■ 



HsvuOt + 2K A sm49 T = 



The boundaries between the various phases in the limit 
C — ► 00 can be found quite simply. The Ila-IIb boundary 
is the locus of points where Ht — -ff/v- From (|l7| ) and 
© we get K = (1+^) 1/2 « 1.13. The Ilb-III boundary 
occurs when Hn — 0, i.e., JC = 1. The III-IV phase 
boundary is the locus of points where Hm = —Ht- This 
^) 1 / 2 «0.85. 



(15) gives IC = (1 



which determines 9t- To find 9s, we need only note that 
the 9t = solution to ( |l5|) is valid for large values of the 
external field. We therefore substitute this value into (|^) 
to find 



1 



(fC 2 



1) sin 2 29 



s 



Hsin 2 U s . 



(16) 



The identification 



= 8K 4 (K, 2 - 1). 



follows immediately since, by definition, 9s is very small 
near nucleation. Substitution of ( |l7| ) into (16) gives 



H = Hff cos 2 9s cos 2 U 



(18) 



which is valid so long as 8t — and H < \H 



N 



The case H™ > is relevant to Phases Ha and lib 
where is distinct from H™ . In particular, the step 
angle 9s increases smoothly as H decreases until the lat- 
ter reaches 



Hg° = 



(19) 



when a magnetization jump occurs because (|18J) has no 
solutions for H < 0. The spin configuration just before 
the jump is precisely that of a 90° domain wall because 
9s = 7T /2 and 9t = 0. As noted in Section II, the jump 
occurs because the domain wall depins from the step and 
sweeps across the terrace so that final state has 9{x) — 



D. Other Analytic Results 

This section presents three analytic results that pertain 
to interior portions of the phase diagram. The first is an 
implicit expression for the nucleation field at any point 
in the phase diagram. The second is an exact expression 
for the entire boundary between Phase He and Phase III. 
The third is the leading correction to the Phase II jump 



(17) field He? when the terrace length is finite. 



For the nucleation field, our interest is the first devi- 
ation of the spin configuration from 9(x) = 0. We thus 
expand (|^) to first order in 9: 



J 



dx 2 



(H + 8K 4 )9 - S(x)2K 2 



(21) 



Without the delta function, the appropriate solution to 
©is 



9 = Acosh(y/(H + 8K A )/Jx) 



(22) 



where A is a constant. Similarly linearizing the boundary 
condition ([5]) gives 



2/± 

dx 



L/2 



-2./i 
dx 



2K 2 8 S . 



(23) 



-L/2 



Combining these results yields the implicit formula 



6 



2K 2 + 2y/j{H N + 8K 4 ) tanh 



H 



N 



8K 4 



J 



= 0. 



(24) 



for the nucleation field Hjy. We obtain a more compact 
form by defining a shifted and scaled nucleation field 
from 

H N {K 2l K 4 , J, L) = -8K 4 + 8K 4 H N (IC, C). (25) 
and substituting ( p5| ) into (pi]). The final result 

K = H 1 / 2 tanh (jCH^ ■ (26) 

gives the nucleation field at any point in the phase dia- 
gram. Note the limiting forms Hn — IC/C for ICC — > 
and Hn = IC 2 for ICC — * 00. These are the Stoner- 
Wohlfarth and single step results obtained earlier. The 
line ICC = 1 can be regarded as a crossover between the 
two. We return to this point in Section IV. 

The IIc-III phase boundary is defined by if at = 0, i.e., 
if TV = 1. Substitution of this into (|26|) gives 



K = tanh(£) 



(27) 



which is the equation of the phase boundary drawn in 
Figure ||. 

We turn finally to a calculation of the jump field Hs in 
Phase II for large but finite terrace lengths. In this limit, 
the domain wall dcpinning picture of the jump is appro- 
priate. The calculation is analogous to the computation 
in Section III.C except that the single step formula (|l5| ) 
is replaced by a more general relation between 9s and 9t 
obtained from a variational form for the spin configura- 
tion near Hs- 

Just below (|l~9|), we observed that the single-step spin 
configuration 0{x) just before the magnetization jump at 



ifg° takes the form of a 90° domain wall. That is, 



tan 9 



(28) 



where A = \j8K 4 j 3 . Since 6s — vr/2 at every step, 
an appropriate trial function for a multi-step system is 
obtained by adding together the ± wall configurations 
from (f28) in the form 



tan 9 = tan 9t cosh Ax, 



which becomes 



tan 9$ — tan 9t cosh C 



(29) 



(30) 



at each step. Expanding (|30|) for large C and small e 
ir/2 — 9 s and 9t gives 



(31) 



Performing a similar expansion on (^|) and retaining 
terms to lowest order in if only yields 



H = W H n^ 2C + ±K 4 9 2 T 



when ( pip is used. The jump field 

H s = 2^8K 4 H™e~ c 



(32) 



(33) 



is the smallest value of H for which solutions to ( 32 ) exist 
for some value of 9t- 



E. Numerical Results 

Numerical methods were used to study three aspects 
of this problem: (i) calculation of the hysteresis loops; 
(ii) determination of the C dependence of the jump fields 
for representative values of IC; and (iii) determination of 
the phase boundaries in the phase diagram. 

The hysteresis loops in Figure ^| were computed di- 
rectly from ([j]). For each choice of control parameters, 
the evolution of the stable energy minimum was followed 
by a combination of conjugate gradient (CG) minimiza- 
tion and spin relaxation dynamics. The initial state was 
chosen as the saturated state and the external field was 
reversed in small steps from a large positive value to a 
large negative value. The CG method reliably follows the 
adiabatic minimum until a magnetization jump occurs. 
But when a jump connects local energy minima that are 
far separated in configuration space, the CG scheme of- 
ten predicts an obviously incorrect final state. To correct 
this, CG was used consistently except in the immediate 
vicinity of a jump. When it predicted a jump, the simu- 
lation was backed up and spin relaxation dynamics used 
to find the correct final state. 

The nucleation field is found readily numerically from 
the general formula (p4|). A more elaborate procedure is 
needed to find the jump fields. Jumps in magnetization 
correspond to discontinuous changes in the spin configu- 
ration. In particular, d9x/dH diverges at both Hs and 
Ht- But since £ is a constant for a given physical situa- 
tion, it must be the case that 



dC 
dH 



dC 
dH 



dC d9 T 
d9 T dH 



= 0. 



(34) 



In this equation, C is regarded as a function of 9t and H 
only since 9s is a function of 9t and H from (||) . We con- 
clude that there is a one-to-one correspondence between 
the divergences of 09T/dH and the zeros of dC/d9x- 

The argument above directs us to find C(9t) for any 
desired choice of IC and H. Once this choice is made, we 
sample many values of 9t in the interval < 9t < 7r/2. 
For each 9t, we solve ([|) for 9s and integrate (0) to get 
C. Figure | shows C{9 T ) for /C=1.25 and H=2, 3, 4, 6, 9. 
The value of H decreases monotonically as the sequence 
of curves is traversed from bottom to top. All the curves 
approach either C — 00 as 9t — ► or possess a semi- 
infinite vertical segment at 9t = that begins at the 
point where the curve hits the left C axis. 
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FIG. 4. Scaled terrace width £ as a function of terrace 
spin angle Qt for /C= 1.25 and H= 2,3,4,6,9. The value 
of Ji decreases monotonically as the sequence of curves is 
traversed from bottom to top. The horizontal dashed lines 
labeled L* ,£s, and Ct are discussed in the text. 
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FIG. 5. Characteristic fields for K, = 1.25. The vertical 
dashed line £ = C* is the I-IIa phase boundary. See text for 
discussion. 



We now argue that the horizontal dashed line labeled 
Cs that is tangent to the local minimum of one of the 
displayed curves defines the physical terrace width for 
which the corresponding value of H is exactly Hs- Hs 
is encountered by reducing the field from large positive 
values where the spin configuration is saturated. The in- 
tersection of the line Cs with the vertical portion of the 
curves for large TC confirms that 8t = at saturation. As 
Ti. decreases, the corresponding curves eventually inter- 
sect the line Cs at small non-zero values of 8t- Finally, 
the intersection occurs at the local minimum of one of 
the curves. This is the curve of Hs because any further 
reduction in field leads to a discontinuous change in 8t 
to the only remaining intersection point on the rightmost 
segment of the C(9t) curves. 

The horizontal dashed line labeled Ct that is tangent 
to the local maximum of one of the curves defines the 
physical terrace width for which the corresponding value 
of H is exactly Ht- But since Fig. [|is drawn for H > 
only, the jump at Ht is encountered by increasing the 
external field from H = where M — 26 1. The inter- 
section of the line Ct with the lowest field curve shown 
confirms that Ot — tt/2. As 7i increases, the curves de- 
velop a local maximum and the intersection eventually 
occurs at this point. This is the curve of Ht because any 
further increase in field leads to a discontinuous change 
in Qt to the only remaining intersection point on the left- 
most segment of the C(9t) curves. 

The evolution of the nucleation and jumps fields as a 
function of C found as described above is illustrated in 
Figure || for JC — 1.25. Figure ||(a) confirms the expo- 
nential dependence of Hs on C predicted in ( |33| ) . Figure 
^(b) shows that Ht ~ C x for the last decade of data 
shown where x — 3-7. 

The relative values of H^, Ht, Hs were used to con- 



struct all the phase boundaries shown in Figure g. Figure 
H is germane to the I-IIa phase boundary. No jump fields 
exist for C < C* and H N > H T > H s for C > C* . 
This is the same terrace length shown in Figure |4| where 
the dashed line C = C* intersects the curve of C(6t) 
for which the extrema (and hence the jump fields) first 
disappear. The I-IIa phase boundary is asymptotically 
vertical as JC — > oo. The limiting value of C*^ is found 
from the same procedure as above by putting 8s = tt/2 
in (0). The result is C*^ w 2.2072. 

Figure fj] shows the C dependence of the nucleation 
and jump fields for JC — 0.5. The absence of the jump 
fields defines the range of Phase I as before. The other 
phases exhibit the relative orderings of the characteris- 
tic fields discussed in Section II, i.e.,. Hn = defines 
the IIc-III boundary and Ht = —Hn defines the III-IV 
boundary. Figure |8| shows the nucleation and jump fields 
for JC — 1.17. The re-entrant behavior Ila — > lib — > Ila 
described in Section II arises because the curves of Ht 
and Hn intersect twice. The transition from Ila to lib at 
fixed C is readily understood. Ht is nearly independent 
of K, because it is related to terrace spin behavior far 
from the steps. But Hm decreases rapidly as K, decreases 
because the torque on step spins is reduced. Eventually, 
Hn drops below Ht for all values of C. We omit a fig- 
ure that shows the Ilc-IIb phase boundary (Hs = H^) 
explicitly. 

We note finally that there is a critical point in the 
phase diagram (fCc,Cc) where Ht, Hs, and Hn are 
coincident. This is the point in Figure || where the I-IIa, 
Ila-IIb, Ilb-IIc, and IIc-I phase boundaries all meet. Our 
best estimate is Kc ~ 1.10 and Cc ~ 0.56. 
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FIG. 6. Asymptotic behavior of the fields Tls and 'Ht — 'H't 
for large C. Note that ffg° = 0.(a) Log-linear plot, (b) 
Log-Log plot. Straight line has a slope of -3.7 
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FIG. 7. Characteristic fields for /C = 0.5. Vertical dashed 
lines denote phase boundaries. See text for discussion. 
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FIG. 8. Characteristic fields for K. = 1.17. Vertical dashed 
lines denote phase boundaries. See text for discussion. 

IV. DISCUSSION 
A. The Reversal Mechanism 

An important conclusion from our analysis is that a 
distinct hysteresis loop topology does not imply a dis- 
tinct mechanism of magnetization reversal. This is im- 
mediately clear from Figure |^ where all four phases are 
present in the CK. -C 1 limit of nearly coherent rota- 
tion and three of the four phases are present in the limit 
of widely separated steps where reversal occurs by do- 
main wall depinning. No sharp transition separates these 
cases. Instead the reversal mechanism smoothly crosses 
over from coherent rotation to domain wall depinning as 
the terrace length or step anisotropy is increased. 

The crossover is most easily understood for the case 
of nucleation which, as noted, always occurs at the steps 
due to the torque exerted on the saturated state by the 
local two-fold anisotropy. When CK ^> 1, nucleation re- 
sults in the formation of a domain of rotated spins around 
each step separated from the unrotated terrace spins by 
a domain wall. Now suppose that C is reduced, say, by 
increasing the vicinality of the substrate. The spins on 
the terrace rotate away from saturation when the domain 
walls begin to overlap. In the limit when C -C K^ 1 , their 
rotation becomes indistinguishable from the rotation of 
the step spins and the coherent rotation picture is a good 
approximation to nucleation. Alternately, suppose that 
K is reduced, say, by increasing the film thickness or by 
adsorbing foreign gases onto the steps. This reduces the 
torque on the step spins so that their angular deviation 
from the terrace spins is not as great. In the limit when 
K <C C^ 1 , this difference nearly disappears and the co- 
herent rotation picture is again appropriate. 

We turn next the first jump field Us- Coburn and 
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co-workers 27 have presented a model of reversal for ul- 
tratliin magnetic films with in-plane magnetization and 
four-fold anisotropy. They assume that the film is well 
described by a single homogeneous domain before and af- 
ter every jump in the hysteresis curve. Domain walls are 
presumed to nucleate at widely separated surface steps 
or other defects. Magnetization jumps occur when the 
energy density gain to make the transition AE is equal 
to a phenomenological energy density e needed to depin 
the wall from the most effective pin in the film. 

This description approximately reproduces our results 
when C is large if we take account of the inhomogeneous 
spin configuration induced by the steps. In Section II, 
Hs was defined as the field when the energy barrier A ow 
vanished. Here, Adw — e — AE where 



AE ■ 



-H/a + 2<r/L+ {A/L)e 



(35) 



and A is a constant with dimensions of energy. The first 
term is the Zeeman energy gain of the saturated state 
compared to the 90° state. The second term is the en- 
ergy cost of the domain walls near the two steps that 
bound a terrace. The last term represents an effective re- 
pulsive interaction between neighboring walls that arises 
from the overlap of domain walls. The terrace spins in 
the overlap region pay anisotropy energy, and the energy 
of the initial state rises compared to the single step case. 
The exponential dependence on wall separation is famil- 
iar from other problems where periodic domains form, 
e.g, the commensurate-incommensurate transition. [ p8| 
The condition Adw — yields the estimate 



H s ~ 2a /L — e + (A/L)e 



-c 



(36) 



This agrees with ([53|) up to the prefactor of the exponen- 
tial if e = 2<j/L. This is not unreasonable because the 
barrier for the two domain walls to depin, sweep across 
their common terrace, and annihilate is associated with 
a spin configuration where the two walls are separated 
by a distance small compared to L but large compared 
to the exchange length W . Of course, e is not distributed 
across the terrace in any physical sense. It is associated 
solely with the particular spin configuration described 
just above. 



B. Comparison to Experiment 

We remarked in Section II. C that the shape of the 
SMOKE loops obtained by Chen and Erskine [[l9) for 
flat and vicinal ultrathin Fe/W(001) appear (to the eye) 
to be very similar to our Phase III and phase II topolo- 
gies, respectively. To see that this is not unreasonable, 
we combine the 25A terrace widths reported in Ref. fl9f | 
with typical values of the magnetic parameters J ~ 10 1 
J, K 2 ~ 1 mJ/m 2 , and K4 ~ 10~ 2 mJ/m 2 || to discover 
that this experiment corresponds to C *~ 1 and }C ~ 1. 
This is indeed in the vicinity of the II-III phase boundary. 



We assigned the same transition to the data of 
Kawakami et al. J2l| for 25 ML of Fe on a sequence of sur- 
faces vicinal to Ag(001). This is still nominally an ultra- 
thin film because the exchange length W = \J J/2K4 ~ 
20 ML using the values above. In fact, the results of this 
experiment lay even closer to the lower left corner of our 
phase diagram than the Chen and Erskine experiment 
because the vicinality is greater. 

The authors of Ref. [£l| analyzed their data with a 
single domain model similar to that of Coburn and co- 
workers p7| except that the step anisotropy was dis- 
tributed across the terraces and the depinning energy e 
was set to zero. Such a model actually yields no hysteresis 
at all-just a magnetization curve with two symmetrical 
jumps. Magnetic parameters were extracted from the ex- 
periment by matching this jump to the average of what 
we call Hs and Ht- In our opinion, formulae similar 
to our ( |l3| ) and (14) for H s and H T should be used to 
analyze the large vicinality data fl29] of Ref. . 



C. Extensions of the Model 

It is easy to think of extensions of the model studied 
here that would render the results more directly com- 
parable to experiment. Probably the most stringent as- 
sumption we make is that the magnetic film smoothly 
coats the vicinal substrate. For relatively small terrace 
lengths, this is possible if the deposition is performed at 
high temperature so that nucleation of islands on the ter- 
races is suppressed and growth occurs in so-called "step- 
flow" mode. (2^] Otherwise, it is necessary to take ac- 
count of the effect of these islands on the hysteresis. This 
was the subject of a previous paper by us |3| for a square 
island geometry and we can use those results to suggest 
the effect in the present case. 

For fixed deposition conditions, island nucleation is in- 
creasingly probable as the terrace length increases. |||] 
For this reason, we focus on the right hand side of the 
phase diagram. The magnetization jump at Hs will be 
interrupted because the domain walls depinned from the 
vicinal steps will not sweep completely across the ter- 
races. Instead, they will be repinned by the channels be- 
tween islands. This introduces additional jump structure 
into the hysteresis curves and likely will alter the coer- 
cive field significantly. We expect little change in Ht 
but there will be an extra magnetization jump before fi- 
nal reversal associated with spins that remain pinned in 
the original saturation direction at island edges perpen- 
dicular to the vicinal step edges. 

The one-dimensional character of our model arises be- 
cause we assumed perfectly straight steps. This is not 
generally the case because the desired step-flow growth 
mode itself induces a step- wandering instability. j3(J This 
instability will have the effect of introducing two-fold 
anisotropics in a variety of directions and a random 
anisotropy model (with spatially correlated randomness) 
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might be a suitable starting point in the limit of large 
waviness. 

Non-uniform terrace widths are another feature of real 
vicinal surfaces that might also be treated in a more com- 
plete model. The result is easy to guess in the pinned 
limit where every terrace acts independently. Otherwise, 
nucleation and subsequent jumps will occur first in re- 
gions of the film with largest step density and eventually 
spread to regions of low step density. 

Finally, we have ignored both perpendicular variations 
in the magnetization and all explicit magnetostatic ef- 
fects. For a vicinal surface, dipole-dipole coupling actu- 
ally induces the spins to lay in the average surface plane 
of the entire crystal |Q rather than in the plane of the 
terraces as we have assumed. When combined with crys- 
tallographic surface anisotropy, this effect induces a two- 
fold anisotropy parallel to the steps at all terrace sites. 
pl[ Such a term is easily included in our basic energy 
expression (Q) and does not appreciably complicate the 
analysis. 

V. SUMMARY 

This work was motivated by the increasing aware- 
ness that the step structure of ultrathin magnetic films 
can have a profound effect on magnetic reversal and 
hysteresis. Our theoretical study focused on perhaps 
the simplest case: a film deposited on a vicinal sur- 
face comprised of uniform length terraces separated by 
monoatomic steps. The magnetization was assumed to 
lay in the plane parallel to the terraces and to vary neg- 
ligibly in the direction perpendicular to the terraces and 
parallel to the steps. We assumed the presence of an 
intrinsic four-fold in-plane anisotropy at every site and 
a two-fold anisotropy at step sites only. Explicit mag- 
netostatics was ignored. Attention was directed to the 
interesting case where one orients an external field per- 
pendicular to the direction of the two- fold axes. The final 
model studied was a one-dimensional, ferromagnetic spin 
chain in an external field with spatially inhomogeneous 
anisotropy. 

The analysis was performed in the continuum (micro- 
magnetic) limit where the spin configuration is repre- 
sented by a function 9{x) that encodes the angular devi- 
ation of the magnetization from the external field direc- 
tion. Four characteristic hysteresis loop topologies were 
found and designated as "phases" in a two-dimensional 
diagram labeled by the natural control parameters of 
the model: a scaled terrace length C and a scaled step 
anisotropy strength K,. 

The hysteresis loops were characterized by a nucleation 
field Hn, where the magnetization first deviates from 
saturation, a step jump field H$ where a jump in mag- 
netization occurs from near saturation to a state where 
many spins are aligned parallel to the steps, and a ter- 
race field Ht where a jump in magnetization occurs to 



the nearly reversed state. For large values of L we found 
H s ~ exp(-£) and H T ~ £~ x with X ^ 3.7. 

In all cases, reversal initiates at the steps because the 
torque applied by the local anisotropy is maximal there 
in the saturated state. No sharp transition separates 
the cases of subsequent spin rotation by nearly coher- 
ent rotation and subsequent spin rotation by depinning 
of a 90° domain wall from the steps. It is a crossover 
phenomenon. The coherent rotation model of Stoner & 
Wohlfarth (SW) is most appropriate in the lower left cor- 
ner of our phase diagram. The step depinning picture is 
most appropriate in the upper right corner of the dia- 
gram. 

To our knowledge, all existing measurements of the 
magnetic properties of ultrathin films on vicinal surfaces 
have been confined to a relatively small portion of our 
phase diagram. We encourage experiments designed to 
explore the remaining terra incognita. 

VI. ACKNOWLEDGMENTS 

The authors acknowledge an intellectual debt to Tony 
Arrott for the the basic premises of this study and they 
are grateful to Lei-Han Tang for his contribution to our 
understanding of the single step limit. Ross Hyman was 
supported by National Science Foundation Grant No. 
DMR-9531115. 



[1] U. Gradmann, in Handbook of Magnetic Materials, Vol- 
ume 7, edited by K.H.J. Buschow (Elsevier, 1993), Chap- 
ter 1. 

[2] B. Heinrich and J.F. Cochran, Adv. Phys. 42, 523 (1993). 

[3] L. Neel, J. Phys. Radium 15, 225 (1954). 

[4] D.S. Chuang, C.A. Ballentine, and R.C. O'Handley, 
Phys. Rev. B 49, 15084 (1994). 

[5] See the various contributions to Ultrathin Magnetic 
Structures, edited by J.A.C. Bland and B. Heinrich 
(Springer- Verlag, Berlin, 1994). 

[6] B. Hillebrands, P. Baumgart, and G. Giintherodt, Phys. 
Rev. B 36, 2450 (1987). 

[7] P. Bruno, G. Bayreuther, P. Beauvillain, C. Chappert, G. 
Lugert, D. Renard, J. P. Renard, and J. Seiden, J. Appl. 
Phys. 68, 5759 (1990); M.T. Kief and W.F. Egelhoff, Jr., 
J. Appl. Phys. 73, 6195 (1993); Y.L. He and G.C. Wang, 
J. Appl. Phys. 76, 6446 (1994); X. Meng, X. Bian, R. Ad- 
bouche, W.B. Muir, J.O. Stom-Olsen, Z. Altounian, and 
M. Sutton, J. Appl. Phys. 76, 7084 (1994); F. Baudelet, 
M.-T. Lin, W. Kuch, K. Meinel, B. Choi, CM. Schnei- 
der, and J. Kirschner, Phys. Rev. B51, 12563 (1995); M. 
Giesen-Siebert, F. Schmitz, and H. Ibach, Surf. Sci. 336, 
269 (1995). 

[8] A. Moschel, R.A. Hyman, A. Zangwill and M.D. Stiles, 
Phys. Rev. Lett. 77, 3653 (1996); R.A. Hyman, M.D. 



11 



Stiles, L.-H. Tang, and A. Zangwill, J. Appl. Phys. 81, 
3911 (1997). 

[9] M.E. Buckley, F.O. Schumann, and J.A.C. Bland, Phys. 

Rev. B 52, 6596 (1995); M.E. Buckley, F.O. Schumann, 

and J.A.C. Bland, J. Phys. Cond. Matt. 8, L147 (1996). 
[10] M. Kolesik, M.A. Novotny, and P.A. Rikvold, Phys. Rev. 

B 56, 11791 (1997). 
[11] E.C. Stoner and E.P. Wohlfarth, Phil. Trans. Roy. Soc. 

A240, 74 (1948). 
[12] B.N. Filipov, Phys. Met. Metall. 22, No. 3, 22 (1966). 
[13] W.F. Brown, Jr., J. Mag. Mag. Mat. 3, 295 (1976). 
[14] A.I. Mitsek and S.S. Semyannikov, Sov. Phys. Sol. State 

11, 899 (1969). 

[15] R. Friedberg and D.I. Paul, Phys. Rev. Lett. 34, 1234 
(1975). 

[16] A.S. Arrott and B. Heinrich, J. Mag. Mag. Mat. 93, 571 

(1991) ; A.S. Arrott, J. Appl. Phys. 69, 5212 (1991); A.S. 
Arrott, T.L. Templeton, and Y. Yoshida, IEEE Trans. 
Mag. 29, 2622 (1993); A.S. Arrott, in Nanomagnetism, 
edited by A. Hernando (Kluwer, Dordrecht, 1993), pp. 
73-85. 

[17] B. Heinrich, S.T. Purcell, J.R. Dutcher, K.B. Urquhart, 
J.F. Cochran, and A.S. Arrott, Phys. Rev. B38, 12879 
(1988). 

[18] A. Berger, U. Linke, and H.P. Oepen, Phys. Rev. Lett. 

68, 839 (1992); W. Wulfhekel, S. Knappmann, and H.P. 

Oepen, J. Mag. Mag. Mat. 163, 267 (1996). 
[19] J. Chen and J. L. Erskine, Phys. Rev. Lett. 68, 1212 

(1992) ; J.-S. Suen and J.L. Erskine, Phys. Rev. Lett. 78, 
3567 (1997). 

[20] W. Weber, C. H. Back, and R. Allenspach, Phys. Rev. B 
52, R14400 (1995); W. Weber, C. H. Back, A. Bischof, 
C. Wiirsch, and R. Allenspach, Phys. Rev. Lett. 76, 1940 
(1996). 

[21] R.K. Kawakami, E.J. Escoria-Aparicio, and Z.Q. Qiu, 
Phys. Rev. Lett. 77, 2570 (1996). 

[22] I.V. Markov, Crystal Growth for Beginners - Fundamen- 
tals of Nucleation, Crystal Growth & Epitaxy (World Sci- 
entific, 1995). The experimental distinction between is- 
land growth and step flow growth for Fe/W(110) can be 
seen easily in the scanning tunneling microscopy results 
of H.J. Elmers, J. Hauschild, H. Hoche, U. Gradmann, 
H. Bethge, D. Heuer, and U. Kohler, Phys. Rev. Lett. 
73, 898 (1994). 

[23] The case of K4 favoring orientations rotated by 7r/4 is 
qualitatively different and will be discussed elsewhere. 

[24] Arrott jl6| argues persuasively that magnetostatics does 
not influence the hysteresis of ultrathin films with in- 
plane magnetization. The present authors have confirmed 
this explicitly for a multilayer version of the island model 
studied in Ref. JH1 using the computational methodology 
of Mansuripur |3l|] . 

[25] This is the equation of motion for a kicked rotor in an 
angle-dependent potential if the variable x is replaced by 
time. This mechanical analogy can be helpful in much of 
the qualitative analysis to follow. 

[26] This is clear for all the hysteresis loops depicted in Figure 
|^ where Ht is labeled except Phase III. There, there is 
an inner loop (not shown) that passes through M = at 
H = 0. 



[27] R.P Coburn, S.J. Gray and J.A.C. Bland, Phys. Rev. 
Lett. 79, 4018 (1997) and references therein. 

[28] P. Bak, Rep. Prog. Phys. 45, 587 (1982). 

[29] Our formulae (|l^) and ( |l4| ) should be modified to reflect 
the actual anisotropies assumed in Ref. [^l[ . Our point is 
that a pure Stoner- Wohlfarth type analysis will be suffi- 
cient. 

[30] G.S. Bales and A. Zangwill, Phys. Rev. B 41, 5500 
(1990). 

[31] M. Mansuripur, The Physical Principles of Magneto- 
optical Recording (Cambridge University Press, Cam- 
bridge, 1995), Section 13.2. 



12 



